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We consider layered superconductors with a flux lattice perpendicular to the layers and random 
columnar defects parallel to the magnetic field B. We show that the decoupling transition tempera- 
ture Td, at which the Josephson coupling vanishes, is enhanced by columnar defects by an amount 
5T d /T d ~ B 2 . Decoupling by increasing field can be followed by a reentrant recoupling transition 
for strong disorder. We also consider a commensurate component of the columnar density and show 
that its pinning potential is renormalized to zero above a critical long wavelength disorder. This 
decommnesuration transition may account for a recently observed kink in the melting line. 



The phase diagram of layered superconductors in a 
magnetic field B perpendicular to the layers is of con- 
siderable interest in view of recent experiments on high 
temperature superconductors [1-4]. Columnar defects 
(CD) induced by heavy ion bombardment provide an 
additional interesting probe [1]. In particular the irre- 
versibility line at low temperatures is enhanced [5,6] while 
within the liquid phase an onset of enhanced z axis corre- 
lation [7-9] was observed. Recent data [10] indicates that 
CD produce a porous vortex matter in which ordered 
vortex crystallites are embedded in the 'pores' of a rigid 
matrix of vortices pinned on the CD's. A sharp kink in 
the melting curve signals an abrupt change from melting 
enhanced by the matrix at high fields to a more weakly 
enhanced melting at lower fields. Theoretical studies on 
CD's have shown a "localization" transition within a bose 
glass phase [11,12]. Recent simulations [13] have been in- 
terpreted in terms of a Bragg-Bose glass with positional 
order which sets in as field increases. Also a "recoupling" 
crossover transition [14] was studied in the vortex liquid 
phase. 

In the absence of CD theoretical studies have shown 
layer decoupling due to thermal fluctuations [15-17] 
or due to disorder [17,18]. At this phase transition 
the Josephson coupling between layers vanishes at long 
scales, i.e. the critical current perpendicular to the layers 
vanishes and superconducting correlations in the z direc- 
tion (perpendicular to the layers) become short range. 
Decoupling involves in principle also proliferation of point 
defects - vacancies and interstitials (VI) [19]. The flux 
lattice is present even in the decoupled phase with the z 
axis positional correlations maintained by magnetic cou- 
plings. In the case of point disorder this phase would thus 
still exhibit Bragg glass type order without dislocations. 

An increase in the critical current at a "second peak" 
transition has been interpreted as due to an apparent dis- 
continuity in the tilt modulus at decoupling [20] . Plasma 
resonance data [21,22] has shown a significant jump at 
this transition, consistent with the decoupling scenario. 
Whether this transition is driven by decoupling alone 



rather than by a sudden dislocation proliferation [23] re- 
mains to be investigated. 

In the present work we consider the effects of CD 
within the flux lattice phase, neglecting VI (whose role is 
discussed below). We find that the decoupling transition 
temperature Td(B) is enhanced by CD. In particular for 
strong disorder the low field form Td(B) ~ 1/B becomes 
Td{B) ~ B at strong fields, hence decoupling followed 
by a reentrant transition into a coupled state, i.e. recou- 
pling, is possible with increasing field. These predictions 
can test wether the "second peak" transition is of a de- 
coupling type. We also allow for a finite component of 
the CD density which is commensurate with the flux lat- 
tice, a component which usually needs to be specifically 
prepared [24]. We find that long wavelength disorder 
renormalizes the commensurate coupling to zero, i.e. de- 
commensuration, above a critical value of disorder. We 
propose that the matrix component of the porous vor- 
tex matter provides a commensurability potential for the 
embedded crystallites. At high fields this enhances the 
crystallite melting temperature, while below the decom- 
mensuration transition the crystallites decouple from the 
matrix, leading to a weaker enhancement of melting. 

We study the classical partition function of L/d 
Josephson coupled layers where L — > oo is the total 
length in the z direction perpendicular to the layers and 
d is the interlayer spacing. The elastic energy of the 
transverse displacement fields u(q, k) in the absence of 
Josephson coupling can be written as [25,26] 

H ^ = irJ2 I (cmq 2 + cUk)k 2 z )\u a (q,k)\ 2 (1) 

k.a J 1 

where a replica index a — 1,2..., n is needed below for 
the disorder average. The elastic constants are [25,26] 
= r/{Ma 2 \ 2 ab k 2 z )\n{l+a 2 k 2 /^), c 66 = T/(16dog), 
where k is the wavevector in the z direction, k z = 
(2/<2)sin(fc<f/2), r = ^ld/(4n 2 \l b ), X ab is the magnetic 
penetration length parallel to the layers, Oq is the area 
of the flux lattice unit cell and <&o is the flux quantum, 
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i.e. $0 = Ba 2 , . Note that the Josephson coupling in- 
duces an additional term in C44 [26], as also shown below. 
The decoupling transition of the pure system (for weak 
Josephson coupling) is given by [16,17]: 



T 



n 



dk 

k c 44(^) 



r 



(2) 



and our principal aim is to obtain the corresponding Td 
in presence of correlated disorder. 

Consider a distribution of CD whose positions within 
a layer are random and uncorrelated. Each of the CD 
has a radius 60 and their average areal density ucd 
is low, ncobQ -C 1. A flux line has a core of ra- 
dius £0 which usually satisfies £0 < b . Once a flux 
line is partially inside a CD it gains its core energy E c 
per layer. The pinning potential per unit area is then 

U pi „{Y) = {E c /(%) J2iP( r ~ r ») witn tl:ic sum on tne CD 
positions and p(r) is a shape function, e.g. p(r) = 1 for 
r < b and vanishes for r > b - The variance, neglecting 
CD overlaps, is therefore 



U pm (r)U pm (V) « E 2 n CD (b /^d 2 (r - r') . (3) 

The average with respect to a flux density involves an 
additional factor (£o/oo) 4 due to the decomposition of a 
sharply peaked flux into harmonics with reciprocal vec- 
tors Q. The replica average at temperature T is then 
[27] 



ab k 



d 2 q 
(27T) 2 



sq 2 L5 k ,ou a (q, k)u b (-q, -k) 



J rf 2 rcos[Q.«(r)-u^(r))]}/2T 



(4) 



where W — E 2 ncD(bo/ao) 4 and only the shortest most 
relevant [27] Q is retained. The cos term above involves 
vectors Q and u° which in the averages below yield 
(Q • u) 2 = Q 2 (u 2 + uf)/2; ui is the longitudinal displace- 
ment which is reconsidered below, but is neglected for 
now as it has no effect on the decoupling transition. The 
parameter s measures a long wavelength random torque 
coupled to a local bond angle [28] 7 = (d x u y — d y u x )/2 
since for transverse modes (Vu) 2 = A^ 2 \ note that the 
usual long wavelength disorder couples to V • u, hence it 
involves only longitudinal modes. 

The long range Bragg glass properties depend on the 
nonlinear cos term in Eq. (4). If this cos is expanded, it 
yields - T, a bl d 2 ru a (r,k = 0)u b (r,k = 0), i.e. a k = 
quadratic term which has no effect on the decoupling 
transition. It is therefore essential to treat the Bragg 
glass nonlinearities properly. 

We also allow for a commensurate term of the CD den- 
sity of the form 

n com = -y c (2d/Q 2 ) J d 2 r cos[Q • <(r)] . (5) 



This term assumes a predesigned component of the CD 
density, in addition to the random one. 

Consider next the Josephson phase, i.e. the relative 
superconducting phase of two neighboring layers. Each 
flux line can be viewed as a collection of point singular- 
ities, or pancake vortices, positioned one on top of the 
other in consecutive layers. Around each pancake vortex 
the superconducting phase follows the angle a(r) which 
changes by 2n in a complete rotation. The Josephson 
phase involves then a nonsingular component 9 n (r) and 
a singular contribution from pancake vortices. The lat- 
ter are positioned at R; + u™ in the n-th layer and at 
R; + u™ +1 in the (n + l)-th layer, where R; is the undis- 
torted position of the l-th flux line. The total Josephson 
phase is then 



0n(r) 



]T[ a (r-R 
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(6) 



where the expansion is justified in the Bragg glass since 
the correlation length in the z direction is 3> d. Wc 
define (including now the replica index) b a (q, k) = 
—2-Kde lkd ^ 2 u a {<\, k)k z /(qa^) so that the Josephson phase 
is 6^(r) + 6°(r). Fluctuations of the 0„(r) field involve 
the Josephson energy as well as magnetic field terms, 



k,a 

- yj J2 I d 2 rcos[C(r) + 6£(r)] 

n,a J 



k)f 



(7) 



where [29] G f (q,k) = 47rd 3 (A afc 2 + k 2 z )/{ T q 2 ). The full 
Hamiltonian is then H = H e i + Hdis + T~L C om + Wj- 

We proceed to solve this system by the variational 
method allowing for replica symmetry breaking (RSB) 
[27,30]. The form of the variational Hamiltonian Ho 
is obtained by expanding the cos terms and replacing 
yj, y c and W by variational parameters zj, z c and 
&ab(k), respectively. The Josephson term involves a 9 , b 
cross term which is eliminated by a shift 9 a (q, k) = 
<9 a (q, k)-u a (q, k)zj(2irk z /a 2 q) exp^/cd/2) / (Gj 1 +zj / d). 
Hence (repeated indices are summed) 



H ° = ^E/ ^{G-J(q,k)u^k)u b *( q , 



k) 



+ [G] 1 (q,k) + z J /d]\9 a (c 1 ,k)\ 2 }, (8) 
G ab(l' k ) = I c 66<? 2 + ca{q, k)k 2 z + z c ]S ab 



-sL—8 k fi - cr ab (k) . 



(9) 



The effect of the nonsingular 9 a is to shift C4 4 (fc) of Eq. 
(1) into 044(9, k), 
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/ j \ o / / \ i , in . at decoupling must be taken before q — > 0. 

C44^<7, KJ — C^[K) - — — — - — -j — j-tt l^"J 

47r(l + A^g + A ab kf) The variational method minimizes the free energy 

Fa + (H — Hq)o where the free energy Fo and the av- 
where X 2 = ^/(16iT 3 zjd). Note that the limit A c — > oo erage (...)o correspond to Hq. This yields 



(11) 



e ab (k) = (WQ 2 /2d 2 T) [ dz[coskzexp(-B ab (z)/2) - S ab V exp(-B ac (z)/2)\ 

zj = yj eM-(T/2) [ [M(^) 2 ( ^ ) 2 G aa {q, k) + (Gj\q, k) + zj/dy 1 }} (12) 

Jq.k Q a G f (q, k) + zj/d 1 

z c = y c cxv{-(TQ 2 /4) [ G aa (q,k)} (13) 

Jq,k 



where J k = J d 2 qdk/(2n) 3 and B ab (z) is given by 

Bab (z) = TO 2 j ^t[G aa ( q , k) - cos{kz)G ab {q, k)} . 

Since the disorder is z independent the off diagonal terms 
B a ^ b (z) are z independent so that <r a ^ b has only a k = 
component; hence RSB is present only at k = 0. It is 
convenient to define G^ 1 (q,k) = J2 b G~ b (q,k) so l na l 
for k ^ G c (q,k) = G aa (q,k). The RSB solution re- 
duces here to a one step form [27], hence G c (q, k) can be 
written with self energies in the form 

G-\q, k) = c 66 q 2 + c 44 (q, k)k\ + z c + E^l - 5 k ,o) + I(k) 

(14) 



£i + z c = (WQ 4 /I6ird 2 c e6 ) exp(-5+/2) . (15) 

B + = TQ 2 J k G c (q,k) is a Debye Waller factor which 
is dominated by large q, k so that c%{k) can be used to 
obtain 



< 16ttT 
\B+\ < In- 



r ~ (T/8da 2 X 2 ab ) + Si + J(7r/d) + z c 
while the function I(k) satisfies for T <C r 



(16) 



2j = 2/je 



■Li¥+[g,w)+(2 



^) 2 G c (g,fe; ZJ =0)]- 1 ]- 1 



The integral is dominated by k 3> g so that for fe > 1/ A a /j 
the effect of the field via Gf(q,k) is negligible for 
Op/SttA^ <C 1. The q integration is then dominated by 
c° 4 (fc) (while C66 serves as a cutoff) leading to a lnzj, 



zj ~ yjexp{- 



(27rd/ag) 2 



87r7 27rc« 4 (fc) + [S 1 + z c + /(A : )]/fc2 



hence 



' d 2 7 fcC ; 
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In z,/} 



(18) 



where an effective elastic modulus is c e JJ (k) - 
[Si + z c + I(k)]/k 2 . The z c term is obvious here by an 
expansion of H CO m Eq. (5). However the Si +I(k) term 
cannot be derived by an expansion of Jidis and a full RSB 
treatment is required. Si acts then as z c and leads to a 
divergence of c e JJ \k) as k — > as already noted in Ref. 
[27]. Since large k dominates the integral in Eq. (18) we 
use I(k) s=s Si + z c so that 



T d « T°[l + (Si + z c ) 



A ah"0 



rln(a /cf) 



(19) 



where T° = rag ln(a /rf)/(47rA2 fc ) from Eq. (2) is the 
transition temperature in the pure system. Since Q ~ 



/r^ Lidiiaiiiuii njmpiji ah m u ill lihj puic D^DLClll. kill. 

j- . 2 [(c 66 g 2 + z c + Si) _1 l/a and W — a^ 4 we have from Eq. (15) Si 



+ Zc 



-G c (g,fc)]. 



(17) 



a . Thus the change in Td due to columnar defects 



is 5T d /T% 



~ -B , up to InS terms. From Eqs. 



Note that for k -> this yields I(k) - |fc| while I(k) w 
Si + z c for large fc, up to logarithmic terms. A condi- 
tion for melting can be estimated by a Lindemann num- 
ber c L w 0.15 [11] so that (ul{v)/al) = B+/Att 2 w c|. 
Hence r is a measure of the melting temperature and for 
T«T B + is small. We note that it is essential to keep 
the nonsingular phase 9 to obtain the correct structure 
factor G c (q,k) in Eq. (14). 

The decoupling transition is determined by the vanish- 
ing of zj. Eq. (12) can be written in the form 



(15,19) we obtain our first principal result, 



2 41 7 — 73\ ~ 1° ( — > n CDKb \ 2{ J) 

r^as m(ao/a) a 



where E c ~ 0.2r [11]. For strong disorder and strong 
fields the CD can dominate and then T d ~ i? increases 
with i?. This allows a reentrant behavior, i.e. for a fixed 
temperature as B is increased a decoupling occurs at 
T d ~ T$ ~ 1/S and then a recoupling would occur at 
a higher field, assuming this field is still below melting. 
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Next we address the commensurability term, which, 
unlike the Josephson coupling, depends also on the lon- 
gitudinal u/(q, k) component. We therefore add longitu- 
dinal energy terms: first an elastic energy of the form Eq. 
(1) with C66, c° 4 replaced by Cn and c l iA , respectively [26], 
and secondly the usual long wavelength disorder coupled 
to V • u [27] which yields the form of the first term of 
Eq. (4) with s replaced by s l . Since a ao originates from 
the W term in Eq. (4) Si and I(k) are common to both 
longitudinal and transverse parts while the location of 
the one step solution changes by a factor cn/(cn + cqq). 
Since c§§jc\\ = Oo/167rA^ h <C 1 the effect on Si is small, 
yet the structure factor for the longitudinal modes (ana- 
log of Eq. (14)) is significantly modified by the same Si 
and I(k). 

The equation for z c depends also on the k — com- 
ponent of G a a(q,k) which involves the long wavelength 
disorder parameters s ,s l . Using inversion methods for 
Gab [27,30] we obtain our second prinicipal result, 

Zc „ y c ( Zc )l/2 (zc) ,Q 2 /(16^ 6 )+ s 'Q 2 /(16^ 1 ) _ (21) 

Si + z c 

Hence at some critical s , s l (where the powers of z c on 
both sides of (21) equal) the commensurability potential 
is renormalizcd to zero. We note that this renormaliza- 
tion is driven by k = terms, i.e. the same derivation is 
valid for a 2D system with point disorder [31]. 

Long wavelength disorder can generate dislocations 
[28] at s > cIqCLq/Wtt, i.e. below the critical value for 
the vanishing of z c . Furthermore, on very long scales 
dislocations will be induced by short wavelength CD dis- 
order as the system is effectively two-dimensional [28,32]. 
We limit our discussion to a Bragg glass domain which 
ignores these very long scale effects. 

The decoupling description neglects point defects, i.e. 
the nucleation of VI. The latter were studied in the ab- 
sence of Josephson coupling and were shown to be gener- 
ated by point disorder [19] , leading to logarithmically cor- 
related disorder for VI . Disordered CD, however, induce 
only a k — component of disorder which has exponen- 
tially decreasing correlations ~ (q 2 + A ah 2 ) _2 (g 2 + Si) -1 , 
i.e. a straight flux line has exponentially decaying in- 
teractions even though each of its pancake components 
has a logarithmic interaction. This type of disorder can- 
not overcome the logarithmic interaction of pancake vor- 
tices; hence, although for any finite (Gaussian) CD dis- 
order the ground state contains a finite density of VI 
perfectly aligned along z, the defect transition tempera- 
ture at which VI uncorrelated between layers unbind is 
not affected by the CD. The true decoupling, which al- 
lows for both Josephson phase fluctuations and for point 
defects, lies in between the above decoupling and the de- 
fect transition. For not too small Josephson coupling the 
transition is near the decoupling one [29,19] and therefore 
the results above should apply. 



Finally, we address the data [10] on the melting curve 
showing a kink at fields Bk ~> B^. Within the proposed 
porous vortex model [10] we suggest that the "vortex ma- 
trix", pinned by the random CD, forms a commensurate 
potential. The lowest harmonic of this potential which 
couples to the flux periodicity has wavevector Q; [har- 
monics with Q' > Q have a (z c /(z c + Zi)) Q ' /2Q factor 
in Eq. (21), forcing a z c = solution]. Since s is a sec- 
ond order effect, we consider s l = W/d 2 ~ B^B 2 , hence 
decommensuration occurs at B ~ B^, the bare propor- 
tionality constant being, however, too small to account 
for the data. The parameters s, s l are relevant parame- 
ters within RG [28] so that their renormalized values can 
be large. The main result is then that elasticity domi- 
nates at large B while disorder dominates at low B, driv- 
ing z c — > 0, in qualitative agreement with the data. We 
propose then to search for an additional phase transition 
line within the solid phase, corresponding to decommen- 
suration, which meets the melting curve at Bk- 

In conclusion we have shown that columnar defects en- 
hance the decoupling transition so that dTd/T® ~ B 2 . In 
contrast, the melting temperature involves the same ra- 
tio, however within a logarithm [see Eq. (16)]; hence 
at weak disorder the enhancement is also ~ B 2 while at 
strong disorder only a weak h\B effect. The B 2 enhance- 
ment at strong disorder can therefore be useful in iden- 
tifying a decoupling transition. Furthermore, for strong 
CD disorder a possibility of a reentrant transition has 
been found, i.e. with increasing field decoupling is fol- 
lowed by recoupling. We have also studied effects of a 
commensurate CD density and shown that its potential 
vanishes above a critical value of long wavelength disor- 
der. This decommensuration transition may account for 
the unusual kink in the melting curve data. 
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